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Using tools from classical signal processing, we show how to determine the dimensionality of a quantum sys- 
tem as well as the effective size of the environment's memory from observable dynamics in a model-independent 
way. We discuss the dependence on the number of conserved quantities, the relation to ergodicity and prove 
a converse showing that a Hilbert space of dimension D + 2 is sufficient to describe every bounded sequence 
of measurements originating from any D-dimensional linear equations of motion. This is in sharp contrast to 
classical stochastic processes which are subject to more severe restrictions: a simple spectral analysis shows that 
the gap between the required dimensionality of a quantum and a classical description of an observed evolution 
can be arbitrary large. 



In Quantum Information Science the dimension of the ac- 
cessible Hilbert space has the character of a resource — larger 
dimensions mean potentially more powerful protocols. Var- 
ious implementations deal with huge Hilbert spaces corre- 
sponding to ensembles of atoms H] or molecules [2], or con- 
tinuous degrees of freedom leading to an infinite dimensional 
space in the first place. But how many degrees of freedom 
are effectively used? Can we assess the dimension of the 
underlying system from observable data without assuming a 
detailed model description beforehand? These questions, par- 
tially motivated by the need of questioning and pinpointing 
the assumptions of security proofs in quantum key distribu- 
tion Q, were recently addressed in the context of 'non-local' 
quantum correlations. 

In H it was shown that a tripartite system can only yield 
violations of certain Bell inequalities of the order of \fd if 
each subsystem has dimension at least d, and in J5J 01 the 
dimension dependence of correlations has been investigated 
in detail for bipartite systems. As these approaches are based 
on static correlations between several parts, the question has 
been raised [5] whether and how the dimension of a single 
system can be observed. 

Similarly, one might want to have a preferably model- 
independent way of assessing the effective dimensionality of 
the systems environment (quantifying non-Markovianity Q) 
or the number of preserved, 'noiseless' degrees of freedom. 

The present work addresses these question from a dynami- 
cal point of view. Given a discrete time evolution of an expec- 
tation value, we ask what can be inferred about the effective 
dimension of the systems Hilbert space or the environments 
memory. We thereby focus on using as little a priori informa- 
tion as possible. When addressing the systems dimensionality, 
the only assumptions are that the evolution is homogenous in 
time and Markovian in the sense that it is performed on time 
scales large compared to the relevant relaxation times of the 
systems environment. If the latter is not fulfilled we will see 
the environmental memory degrees of freedom in our dimen- 
sion count. 

For this analysis two standard tools from classical signal 
processing [8| can be employed: delayed embeddings and 



analysis in the frequency domain. These will allow us not 
only to tackle the above question but also to address the con- 
verse: which Hilbert space dimension is sufficient for a given 
sequence of measurements? and to compare the efficiencies 
of quantum versus classical descriptions of a given evolution. 
While every sequence produced by quantum mechanical evo- 
lution can in principle be described by a classical stochastic 
process, we will easily see that quantum mechanics can be 
arbitrarily more efficient in terms of the required number of 
discrete degrees of freedom. 

Not surprisingly we will see a connection between the ef- 
fective dimension and the number of conservation laws and 
get a glimpse on more (in particular spectral) information 
which can be obtained from the evolution. 

PRELIMINARIES 

Our interest lies in the discrete time evolution of expecta- 
tion values of the form 

(A(t)) = tr [AT\p)] , (1) 

where p is a density matrix acting on a d-dimensional Hilbert 
space, T is quantum channel, i.e., a trace-preserving com- 
pletely positive linear map with equal input and output space 
and A = A^ an observable. For simplicity we will in 
the first part assume that ([TJ is a half-infinite sequence, i.e., 
t £ No = {0, 1, . . .}; the extension to finite sequences and 
noisy data will then be discussed at the end. 

Note that the description in ([T} assumes homogeneity in 
time and Markovness in the sense that future evolution only 
depends on the state at present and not on its history. This 
means that given the l.h.s. of ([T} the Hilbert space underly- 
ing the description of the r.h.s. has to contain all effectively 
relevant degrees of freedom. Hence, if a system of dimension 
ds undergoes a non-Markovian evolution due to memory 
degrees of freedom in the environment, then d = d$ + dE- 

It will sometimes be advantageous to consider p and A as 
elements of a d 2 -dimensional Hilbert space H equipped with 
the Hilbert-Schmidt scalar product (A\p) := tr As a 
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linear map T has a matrix representation on TL which we will 
denote by T, Using matrix units as a basis of TL we can write 
T = J2k Kk^Kk, where {Kk} is a set of Kraus operators of 
T i) = J2k K k ■ K\. While T refers to time evolution in the 
Schrodinger picture, we will denote by T* the respective map 
in the Heisenberg picture so that tr [AT{p)\ = tr [T* (A) p] . 
We will denote by 

C := {H = H^\(H(t)) independent of t G N } (2) 

the space of conserved quantities which obviously includes all 
H = T*(H) and in particular 1 € C as the evolution is trace 
preserving. Note that C depends on T and p. 

A quantum channel T will be called ergodic w.r.t. a state p 
(an observable A) if the orbit generated by T* (T**) spans the 
entire space of d x d matrices. 

ASSESSING THE DIMENSION 

The central tool in this section is the space V spanned by all 
delayed vectors of the form 

v t =((A(t)),(A(t+1)),...), re N . (3) 

The employed approach, often called method of delays, is par- 
ticularly widespread in the analysis of chaotic dynamics |9| 
and it provides the following simple and tight relation: 

Proposition 1 (A bound on the dimensionality) Consider 
the space V = span{v T } T( zm Q spanned by the delayed vectors 
obtained from a sequence of the form (|7J. Then 

dimC + dimV < d 2 + 1, (4) 

where d is the dimension of the underlying Hilbert space and 
dim C the number of linearly independent conserved quan- 
tities. Equality holds in ^ if T is ergodic w.r.t. A and 
dim V — d 2 iff it is ergodic w.r.t. A and p. 

Proof. Consider a basis {Hi}, i = 1, ...,£) := dimC of the 
space C of conserved quantities. Then {Hi — tr \pHj\ 1} span 
a D — 1 dimensional subspace of TL which is orthogonal to 
the space spanned by T t \p), t S No. Hence the latter space, 
denote it by Sj, has dimension d 2 — D + 1 and 

(A(t)) = {A Si \T%\p Si ), (5) 

where the index 3 refers to the restriction onto Sj. Since the 
minimal polynomial |10| of has degree at most dim.f) 
there are complex coefficients Cj such that 

d 2 -D 

3=0 

Recalling that (v T )k — (A(t + k — 1)) this implies that there 
are at most dim Sj linearly independent vectors in V since 

d 2 -D 

(v d 2_ D+x ) k = 22 Cj (vj) x+ k-i Vx,keN. (7) 
3=0 



Let us now assume that T is ergodic w.r.t. A, which means 
that {(A|T*} te fsj spans Ji. Suppose that Q would not be an 

equality. Then ^2^ =0 D CjVj = for some c which implies, 

by ergodicity, that J2f=o° Cjfi+ n \p) = for all n G N . 
Therefore dim Sj < d 2 — D so that for the orthogonal comple- 
ment dimf)- 1 > D. However, this contradicts the inequality 
dim Sj < D — 1 which comes from the fact that every ele- 
ment in Sj 1 - is in C and in addition 1 e C \ Sj- 1 . So, ergodicity 
w.r.t. A implies equality in Q. 

If T is in addition ergodic w.r.t. p then Sj = Ti implies 
D = 1 so that indeed dimV = d 2 . We finally prove the 
converse again by contradiction: assume linear dependence 

d 2 — 1 

of the form 5Jj=o C 3 V 3 ~ 0- Then for all a, b € No we 
have (A\f a (J2j c j fi)f b \p) = so that, due to ergodicity, 
T would have a minimal polynomial of degree less than d 2 . 
However, this would again imply the existence of a proper 
subspace S) contradicting the assumption dimV — d 2 . To 
see this recall that a minimal polynomial of smaller degree re- 
quires an eigenvalue A with geometric multiplicity larger than 
one ifTOl . Denoting by \(j>) a linear combination of the corre- 
sponding left eigenvectors, we get (<f>\T t \p) = X t (cj)\p) such 
that there is always a <j> J_ Sj. □ 

Some remarks. Depending on the assumptions we may use 
the above result for different purposes: (i) assuming Marko- 
vianity it provides a lower bound for d, (ii) assuming we know 
d in addition it yields an upper bound on dim C and (iii) if we 
only know the dimension of the system d$ = d — dE it gives 
a lower bound on the number dE of effective memory degrees 
of freedom in the environment. In fact, if (A(t)) exhibits alge- 
braic decay then, as one would expect, Eq.Q leads to d = 00. 

Prop. [Tj is easily generalized to the case where, instead of 
a single expectation value, we observe a set of observables 
{A a } or, equivalently, take higher moments of the observable 
into account (i.e., A a — A a ). Then the delayed vectors have 
to be replaced by 'delayed matrices' so that (A a (r)) are the 
entries of the first column of the matrix v T . Eq.Q is then 
obtained in just the same way where dim V is now the number 
if linearly independent matrices. 

QUANTUM EVOLUTION FOR GIVEN SEQUENCES 

The previous section provided a lower bound on the dimen- 
sion of the Hilbert space in terms of the dimension of the space 
V of delayed vectors which are in turn solely based on (A(t)). 
Clearly, there cannot be an upper bound to the 'true' dimen- 
sion since there might always be redundant degrees of free- 
dom. Nevertheless, we can provide a converse to the above 
observation in the sense that there always exists a quantum 
representation in a Hilbert space of dimension not much larger 
than dimV. Remarkably, such a converse does not exist for 
classical evolutions (see subsequent section). 

Proposition 2 (Quantum representation) 

Given any bounded sequence (A(t)) G R, t £ No there is 
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FIG. 1: a) Eigenvalues of the evolution operator lie on the unit disc 
and can be obtained from the observable time dependence of any ex- 
pectation value. For a classical description they have to be located 
in a region depending on the number d c of degrees of freedom. For 
d c = 2 this is the real line, for d c — 3,4 the dark and light grey 
regions become additionally accessible, b) A simple damped oscilla- 
tion leads to eigenvalues e _7± *", implying that a quantum mechan- 
ical description can access the entire unit disc even for d — 2. 



always a quantum state p, a quantum channel T and a Her- 
mitian observable A acting on a Hilbert space of dimension 
dim V + 2 such that ([7]) holds. 

Proof. We begin with the fact proven in Lemma [T] in the 
appendix: there is always a contractive matrix M and vec- 
tors R, L of dimension dim V such that (A(t)) = (L|M'|i?). 
We will proceed in two steps: first establish complete posi- 
tivity by adding one degree of freedom and then impose the 
trace preserving property by adding another degree. Define a 
square 'Kraus-operator' C = 1 © M of dimension dim V + 1 
and, referring to the same block structure (C © C dlmV ), a 
vector |*) = |0 © R) + |1 © 0) and an 'observable' B = 
(|1 © 0)(0 © L\ + h.c.)/2. Then, using that (A(t)) e M we 
have (L|M*|i?) = tr [B C' t |*){*|C7+ t ] . In order to make 
this trace preserving we embed it again, now referring to the 
block structure C © C dimV+1 . With K := © C, A := 



(0 (B B)\\^\\ 2 , p := (0© we obtain that 

T(p) = KpK^ + |1 © 0)(1 © 0|tr [(1 - K*K)p] (8) 

indeed satisfies ([TJ for the chosen p and A, and since 
1 1 M | |oo < 1 implies K^K < 1, T is a valid quantum channel. 

□ 

We do not know whether there are more efficient construc- 
tions. Clearly, from Prop{T]we know that one needs at least 
d > \/dim V, but this will most likely not be achievable in 
general. 



SPECTRAL ANALYSIS AND SEPARATION FROM 
CLASSICAL DESCRIPTIONS 

So far we investigated the dimensionality of the system 
based on the sequence (A(t)). The information obtained can 
be refined when going to the frequency domain by considering 
the function C : C — > C (the 'z-transform' [8 1) defined by the 
series 



C{z) 



1 



(A{t)) 



(9) 



tSN 



This converges outside the unit circle and can be define inside 
by analytic continuation. In this way we obtain for a sequence 
of the form ([TJ 



C(z) =tr A(z id - T) (p) 



(10) 



so that poles of L correspond to eigenvalues of T. Note that 
the latter lie always inside (or on) the unit circle, there is 
an eigenvalue 1, and complex eigenvalues come in conjugate 
pairs. While there are restrictions ifTTl for instance for the de- 
terminant, i.e., the product of eigenvalues, quantum mechan- 
ics does not impose any further constraint on the location of 
eigenvalues: any point on the unit disc is possible even for 
d = 2. The simplest example for this is a damped Rabi os- 
cillation leading to (A(t)) = e~ 7 * cos uit with poles of C at 

It is instructive to compare this with a potential classical 
description of the sequence (A(t)). So assume there are d c 
states to each of which we assign an initial probability pk, 
k = 1, . . . , d c . The evolution of these probabilities for a single 
time-step is governed by a stochastic matrix S and in the end 
a measurement outcome £ R is assigned to the fc'th state. 
In this way we arrive at 



{A(t)) = (a\S t \p). 



(11) 



Yet, the poles of C correspond to eigenvalues of S, which 
share the basic properties mentioned above. However, the 
classical description imposes additional constraints on the lo- 
cation of the eigenvalues depending on the dimension d c . In 
particular, they have to be located inside the convex hull of all 
roots of unity up to order d c . That is, the unit disc will not be 
entirely covered for any finite d c . A more complete character- 
ization of the location of eigenvalues is given in |[T2l [131 and 
shown in Fig{T]for d c = 2, 3, 4. 

A simple consequence of this analysis is that in terms of 
the required degrees of freedom a quantum mechanical de- 
scription of a sequence (A(t)) can be far more efficient than 
a classical one. In the above discussion the separation be- 
tween a quantum and a classical description comes from the 
simple fact that oscillations are easier to describe in terms of 
probability amplitudes than by using probabilities. Certainly 
more sophisticated separations can be found, however, a com- 
plete determination of d and d c from a given sequence (A(t)) 
seems to be a daunting task (despite considerable results on 
the classical side, cf. lfl4l[T5l ). 
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FINITE AND NOISY DATA 

So far we addressed the ideal case of a half-infinite and 
noiseless sequence — also noiseless in the sense that the ex- 
pectation values are known precisely which requires infinite 
statistics. It is, however, straight forward to analyze finite and 
noisy data along the same lines. Let us begin with a finite se- 
quence (A(t)), t < 2(N — 1) and consider the N x N matrix 
Vki '■= (A(k + 1 — 2)). As the rows of V are a truncation of 
the vectors v T , r = 0, . . . , N — 1 we have 

dimV > rank^, (12) 

with equality if N > dim V. If the data are noisy or suffer- 
ing a significant statistical error, then V will typically be of 
full rank. However, if an error estimate is available we may 
consider an effective rank of V by disregarding all singular 
values below a certain noise threshold which is set by the es- 
timated amount of errors. More precisely, assume that V is 
perturbed by some V e (i.e., we actually observe V' = V + V e ) 
where \\V e \\ < e. Then, by application of the singular value 
inequality [16|, we get 

rankV > min{/c | s k+1 (V') < e}, (13) 

where si(V) is the Tth largest singular value of V. 

In Figj2] the behavior of these singular values is de- 
picted graphically for an example of a unitary evolution 
tr \ AU t pU^] with d = 3 and randomly chosen p, A, U. As 
U preserves its eigenstates and has (due to the random choice) 
no other conserved quantities, we have dim V = 7 which is 
well reflected in the singular values of V for small enough er- 
rors. In fact, in such unitary examples the dimension estimates 
appear to be surprisingly stable up to errors which make up a 
considerable fraction of the signal. 

For the spectral analysis finite and noisy data seem to be 
more involved to handle. It is not difficult to show 1 10 1 that (|9]l 
is a rational function with exactly dim V poles and hence we 
are facing the problems of (i) giving a good (the best) rational 
approximation (with a fixed number of poles) of a function 
given by its Taylor coefficients and (ii) obtaining the poles 
of a rational function given by its (noisy) Taylor coefficients. 
The abstract solution to (i) is given by the Adamyan-Arov- 
Krein Theorem [ 17 |. In practice (i.e., for finite data) a possible 
approach is the use of Pade approximants ifTHIl . For (ii), one 
can use for instance the QD-algorithm of Rutishauer lfT9ll . 

DISCUSSION 

In this letter, using tools from classical signal processing, 
we have discussed how to determine the Hilbert space dimen- 
sion needed to explain observed data of an evolving quan- 
tum system. This introduces a new paradigm, beyond Bell 
inequalities and the analysis of correlations, to obtain such es- 
timates with minimal assumptions in the model, in our case 
homogeneity and Markovianity. In particular, the method can 
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FIG. 2: top; a finite sequence (A{t)), t < 100 obtained from a ran- 
domly chosen unitary dynamics of a spin-1 quantum system (d=3). 
As there are 3 conserved quantities (the eigenstates of the evolution) 
we have dim V < 7. bottom: log-plot of the 15 largest singular val- 
ues Sj of the corresponding 50 x 50 matrix V as a function of the 
standard deviation of added Gaussian noise (1%-10% of the signal). 
While for small noise the largest 7 singular values are clearly sepa- 
rated by a threshold, this washes out as the error increases — for too 
much noise the data could as well be explained by smaller dim V as 
expressed quantitatively by Eq.l|13}. 



be used for single systems, which answers a question posed 
in (for a different static approach, based on several in- 
put states, see the recent work |20|). We have also seen that 
we can use the method to quantify the non-Markovianity of 
an evolution via the effective dimension of the environments 
memory — complementing the static approach of [7]. Finally, 
the analysis of spectral information revealed a dramatic differ- 
ence between the dimensions needed to give a quantum resp. 
classical explanation of the data. Clearly, this does not outlaw 
a classical description (as the Garg-Leggett inequalities EH 
do, albeit based on an extra arguable assumption) but it shows 
that quantum mechanics can be much more efficient in terms 
of number of degrees of freedom. For the future, it would be 
nice to combine the two approaches, evolutions and correla- 
tions, to obtain more information. 
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Appendix A 

Lemma 1 (Linear representation) For every bounded se- 
quence (A(t)) G C, ( £ No there are vectors of 
dimension dim V and a respective matrix M with \ \M | < 1 
such that 



(A(t)) = (L\M t \R). 



(14) 



Proof. Consider N x N matrices V M := (A(k + 1-2)} 
and VL :— (A(k + I — 1)) with N large enough so that 
rankF = rankV' = dimV < N. Let V = V L V R be 
a singular value decomposition where Vr,V^ have dimV 
rows. Then \R) := Vr|1), (L\ := (1\Vl and (using pseudo- 
inverses) M := V ] ^ 1 V / V Ii 1 satisfy Eq.(14i obviously for 
t < 2N — 2. For larger values of t, however, the same re- 
lation has to hold by linear dependence, i.e., the fact that by 
construction dimV = rankV. boundedness of the sequence 
(|(j4(f))| < const.) implies that the spectral radius of M is 
at most 1 and that eigenvalues of magnitude one have one- 
dimensional Jordan blocks (i.e. their geometric multiplicity 
equals the algebraic multiplicity). To see this first note that 
powers of a Jordan block J(A) are Toeplitz matrices with 



first row J(X)\ k 



Hence for |A| 



1 



(l\ J(A)*|r) is a polynomial in t which remains only bounded 
for growing t if it is zero. 

Hence M can only have Jordan blocks for |A| < 1. Let 
us finally show that M is therefore similar to a contrac- 
tion: write J(X) = XI + N with N a nilpotent matrix with 
ones on the first upper diagonal. Then J(A) is similar to 
J' := XI + (1 — | A| )iV and using the triangle inequality for 
the norm 1 1 J' 1 1 oo < 1 ■ As the representation of the sequence is 
preserved under similarity transformations of M (and respec- 
tive transformations of R, L) we can always w.l.o.g. chose M 
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